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TECHNICAL MEMORANDUM 


COMPUIERIZED DATA COLLECTION AND REDUCTION 
FROM AN X-RAY DIFFRACIOMEIER 


INTRODUCTION 


A scries of computer pro^^rams has been written for use with a Philips X-ray 
Diffractometer and u Hewlett Packard 9825A Desk Top Computer. These programs 
permit the coiiection and storage on disk of tl\e number of X-ray counts and the 
associated 2 theta angles across line profiles. Automatic background subtraction, 
integrated intensity, correction for the angular dependence of the lorent*, polarisation, 
and atomic scaittering factors, peak location. K alpha 2 removal, and calculation of 
I'uuricr coefficients are performed. This technical note is documentation for these 
pi\)grams aiul should provide a guide to tlieir use. The programs are written in HPL 
which is a Hewlett Packard variation of BASIC. The programs are written for a 
ctjmputor configured with a disk drive, but they can be easily modified to run from 
the cassette tape drive that is integral to the 9825A Computer. 

Copies of these programs arc available on cassette from tlie author. 


HARDWARE 


With the exception of the Angle Mode Programmer (AMP) and tlie Output Inter- 
face. tlu' Philips Electronics ciiuipment used is stmidard X-ray diffraction hardware. 

'1 he opcratioti of the AMP and Output Interface is as follows: data initiation is 

started by pushing "run" on the AMP and the diffractometer advances to the first 
angle that has been programmed into the AMP. This angle is sent to the Output 
Interface. Control is turned over to the Scaler/Timer. After a preset length of 
time the Si'aler /'rimer stops data act|uisltion . sends the number of counts that it has 
collected and the length of time that it collected data to the Output Interface, and 
sends a signal to th.e A^lP that it has completed data acquisition. The AMP causes 
the diffractometer to advance Viy the step size set on the Motor Conti’oller and com- 
pares tlie new angle of the diffractometer with the scan stop angle that has been 
prep rt> gram med into t!ie AMP. If the new angle of the diffractometer exceeds the 
scan stop angle, data acquisition is haltetl, and the measurement is completed. If 
the new angle dues not exceed the scan stop angle, the data atxiuisition cycle is 
repeated . 

During this cycle the Output Interface has received the number of counts 
collected during tlu' preset time, the length of the preset time, and the angle of the 
diffractometer. This data is alivady digitizeil going into the interface. The Output 
Interface appmpriatel y formats tlus string of data for the kind of output device that 
is used. The formatting controls sliould be .set as follows: 


M ode — repeat 
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Parity - - odd 

Format — normal 

Output --- TTY. 

With this formatting a typical output appearing' on ttie DEC writer is shown in 
Figure 1. Column A is the number of counts <»Hectf;<i. Column B is 100 times the 
number of seconds in which the counts were collect e • Coliunn C is 100 times the 
diffractometer 2-theta angle. Column D has no signifiiArce here. 

To interface the Hewlett Packard 9825A Computer to the Philips equipment a 
Hewlett Packard 983036A Option 001 Serial Input/Output (I/O) Interface is required. 

Pin 1 is connected to instrument ground, and pin 3 is connected to the output sign'll 
of the Philips Output Interface. Since the Philips data output is BNC and the Serial 
I/O interface is a 25 pin connector, these connectkms must be made carefully to avoid 
bending or shorting pins. Certain switches internal to the I/O interface must be set 
for proper operation. The switches must be set to: "Current Loop," 2 stop bits. 

Odd Parity, and 8 bits/character. The device select code is 11 and can be set on the 
outside of the I/O interface. 

The basic 9825A computer accepts Read Only Memory (ROM) chips that extend 
the software capabilities. The programs described in this report assume that the 
Advanced Programming, Advanced I/O, String Variable, and Disk Control ROMs are 
in place. The Hewlett Packard 988^ Flexible Disk Drive must also be connected to 
the computer. The 9825A has an integral cassette tape drive that can be used in 
place of the disk drive. During line by line documentation of the programs, some of 
the program modifications that are required to operate this cassette drive will be 
mentioned. Generally, disk operation is preferable because of more flexible addressing, 
greater capacity and speed, and a lower error rate. 


SOFTWARE 


"datrd" 

This program is stored on track 0, file 1 of the cassette and reads pairs of 
X-ray counts and the diffractometer 2 theta angle, prints out this data, and stores 
the counts and angles on disk for future use. A listing of the program is in Figure 
2. The following discussion of the program should help the user to understand how 
the program works and suggests possible modifications that may be useful for par- 
ticular applications. The numbers refer to lines in the program. 

1. The user is asked by the computer how many data points are to be collected. 
The number of points, of course, depends on how wide the line is and the 
step size between angles. 

4. D is a dummy variable and is not used except as a storage for serially 
output data. 

5. C is the vector containing the X-ray counts. 
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7. B is 100 times the number of seconds in which data is collected. If opera- 
tion with variable time increments *s desired, the counts vector can be 
normalized by dividing by B to give counts per second instead of just 
counts. 

8. A is the vector containing 100 times the 2-theta angle of the diffractometer. 

10. The correct 2- theta angle is printed. 

14. The names of the files under which the data are stored are listed. The 
first word after "files" is the name of the file under which the angles are 
stored, and the second word after "files" is the name of the file under 
which counts are stored. These names must be entered before the program 
is run and space must be reserved on the disk for them using an "open" 
command. Any previous information in these files will be lost. 


"Xfoure" 

This program reads data off the disk, performs a background subtraction, 
computes tha integrated intensity of the line, corrects the line for the angle depend- 
ence of the polarization factor and the atomic scattering factor, calcxilates the peak by 
a three point parabolic fit, subtracts the K-alpha 2 component, again calculates the 
peak, and finally calculates the Fourier coe^'ficdents normalized to the area of the peak. 

Except for the Fourier coefficients all these data are standard information that 
the X-ray diffract ionist often wants to know about a given line. The Fourier coeffi- 
cients are used to deconvolute a given profile from the line profile of a well annealed 
specimen that contains no stress and has grains large enough to cause no particle 
size line broadening [1,2]. If the asymmetric coefficients are small compared to the 
symmetric coefficients (as they usually are), the trtie line profile (that is not includ- 
ing instrumental broadening) can be calculated by dividing the observed coefficients 
by the same coefficients of a well annealed specimen. The quotients are the Fourier 
coefficients of the line due only to various X-ray active defects such as inhomogeneous 
stress, small particle size, dislocations, vacancies, and stacking faults. More details 
can be found in Warren [3]. 

"Xfoure" is located on track 0, file 0, and the listing is shown in Figure 3. 
Important lines will be commented on inaividually. 

4 & 5. The computer asks which disk files to use for the angle data and the count 
data. Tliis data must have been previously stored on the disk under the 
appropriate file name. If operation of the tape cassette is desired, lines 4 
through 9 should be replaced by commands to load tine angle data into 
vector A and the counts into vector C. File names cannot be used so the 
program changes must be made before the program is run. 

10 to 16. The data is formatted and output for the operator to verify. 

18. The computer asks the operator to examine the output data and determine 
where the line begins and ends. It assumes that the angles are entered in 
increasing 2 theta values; "L" is the lefthand value where the line begins, 
and "R" is the right hand value as though one were looking at a strip chart. 
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22. A linear fit is made to the background. 

23. Any negative count values resulting from statistical fluctuations in the line 
tails are set equal to zero. This is the only explicit smoothing in the pro- 
gram. A standard multipoint smooth was not used because it reduces 
angular resolution, causes periodicities that may appear in the Fourier 
coefficients, and some smoothing occurs during deconvolution in any case. 

25. If this line is removed, the background subtracted line will be printed out. 

28. The correction for the Lorentz and polarization factors starts here. The 
Lorentz factor is a geometric factor caused by the size of the diffracted 
cone of X-rays, and the polarization factor exists because electrons scatter 
X-rays preferentially in the forward direction so parts of a line profile on 
the low angle side are less than on the high angle side. 

32. This line is only correct for a copper X-ray tube and a graphite mono- 
chromator. If some other conditions are used, ".75 cos Afl]" should be 
replaced by "cos 2 2a cos AH]" where 2a is 180 degrees minus the diffraction 
angle of the monochromator. If no monochromator at all is used this line 
should be replaced by 


(1 + cos A[l] 

Mli A[I] 
2 ® 2 


Q 


35. If this line is removed, the Ixjrentz- polarization factor corrected line will be 
printed out. In this program the LP correction is made before the Fourier 
coefficients are calculated. At large values of theta the main cause of 
broadening in the line of an anneated specimen is spectral breadth ot the 
X-ray source, while at small theta values the main cause of broadening is 
due to instrumental misalignments. If the broadening caused by misalign- 
ments is dominant, the LP correction should be made after the Fourier 
coefficients are calculated [4]. In any case the differences are small. 

40. The Atomic Scattering factor correction subroutine is called. 

47. Q is the corrected sum of all the counts in the line. This is the integrated 
intensity. When comparing the integrated intensities of different lines, the 
angular interval between points must be corrected for. This interval is 
pi'inted out by line 48. 

51. The numbers of three equally spaced points that span the peak are asked 
for and must be entered. 

53 to 55. The corrected number of counts that occur at these points are printed out. 
For best accuracy these points should be 85 percent of the peak above 
background (51. 

59. A parabola is fit through these three {x>ints and the peak of the parabola is 
calculated [6], Kxperience has shown that the peak position is not very 
.sensitive to which three points are chosen. Typically, the peak can be 
repeated to ±0.002 degrees if the 85 percent rule is observed. Naturally, the 
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71. 
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94. 
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statistically uncertainty of the data must be very small for such accuracy. 
This can be assured by long counting times and Is not usually a problem 
near the peak where the count rate should be high. Another source of 
statistical inaccuracy comes from using samples with tco larg^ a grain size 
to give a truly random distribution of particle orientations. Usually the 
particles should be less than 15 ym for cubic crystals and smaller for 
materials with a lower degree of symmetry [7]. 

The K-alpha 2 component of the X-ray line will now be removed. 

rl and r2 are the wavelengths in angstroms of the K-alpha 1 and K-alpha 2 
radiation. The numbers entered in the progiam are for copper radiation. 

If another tube is used the appropriate values must be entered. 

The angular separation of the components of the doublet is calculated at the 
low angle side of the line and for a weighted average of the wavelengths 
from the tube. 


Tliis angular separation is printed out. 

The data occurs in discrete step intervals which are not likely to be even 
submultiples of the angular separation so some interpolation between steps 
will be necessary. This line calculates the whole number of steps in the 
angular separation and outputs this number in line 68. 

The angular separation is calculated for each angle across the line. This 
gives slightly higher accuracy than the standard Rachinger correction 
which usually calculates the separation once for the peak and asstames that 
this separation is constant across the line [8]. To obtain three decimal 
place accuracy from data taken at intervals of 0.01 degrees, this higher 
precision method is necessary. The assumption is still being made, however, 
that the K-alpha 1 and K-alpha 2 peaks have the same shape, which is not 
strictly true. 

Here is the actual Rachinger correction to the number of counts. For inter- 
polation purposes the line is assumed to be linear between data points, and 
the K-alpha 1 line is assumed to be twice the size of the K-alpha 2 line. 

The linearity assumption is quite good for the small step sizes (0.01 degree) 
usually used. The assumption about the ratio of the peaks is also good if 
a low resolution /high reflectivity monochromating crystal such as the Philips 
graphite crystal is used. A sharper monochromating crystal may alter this 
ratio. If it is desired to use a different ratio, the factor ”0.5" which 
occurs twice in this line should be replaced by the correct ratio. 

If this line is removed, the K-alpha 1 line will be printed out. 

The K-alpha 1 peak is printed out using a three point parabolic fit. 

N is the number of Fourier coefficients that will be calculated. Here it is 
set to 12 but can be easily changed. 

The line is moved so that the peak occurs at the origin. This is not 
strictly necessary for the Fourier analysis but it simplifies writing the 
Fourier analysis subroutine. 
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102. The Fourier analysis subroutine requires an odd number of equally spaced 
points in the line. This line determines if there is an odd number of points in the 
line and adds a point with 0 counts on the hig^h angle of the line if needed. 


104. The vector Y now contains the count information for the K-alpha 1 line. 

106. If this line is removed the vector Y is printed out. 

114. The constant term in the Fourier analysis is printed ouc as A (0). 

115 to Formatting for the Fourier coefficients. In line 117 the coefficients are 

118. normalized by the constant term. The constant term, A (0), is the area 

under the curve which is not relevant for line broadening analysis. By 
normalizing the coefficients the shapes of lines collected at different X-ray 
beam currents and for different time intervals can be compared. The sin 
and cos coefficients are output with a prime (') to indicate this normaliza- 
tion. A' is the cosine coefficient and B' is the sine coefficient. X-ray 
lines are nearly symmetric so B' is much smaller than A’ coefficient. The 
B' coefficients are not very accurate because these coefficients are sensitive 
to changes in the peak position as small as 0.002 degree. This is not 
surprising since the line asymmetry is small to start with and any change 
in asymmetry is likely to have a large effect . 

120. The calculation of Fourier coefficients iising Simpson's rule to evaluate the 
integrals begins here. 

123. X and Y are the current value of the angle and the count. 

124. Except at the end points of the evaluation, M alternates between 2 and 4 in 
Simpson's rule. 


127. The values that go into the summation in Simpson's rule are calculated here. 


129 to 
132. 


The constant term A (0) is calculated. 


135. The atomic scattering factor correction is started here. This correction 

arises because atoms scatter X-rays most efficiently in the direction of the 
X ray lx?am. The scattering cross section varies inversely with the sine 
of the scattering angle. Since this factor is angle dependent, the shape 
of the line is changed and must be corrected for. The atomic scattering 
factor depends in a rather complicated way on the wavelength of the radia- 
tion and on the atomic number of the atom tliat is doing the scattering. 
Values of the atomic scattering factor are available and can be used for an 
alloy assuming an average atomic number for the alloy. If the components 
in the alloy are close in atomic number, tliis averaging is a good approxima- 
tion. If the components of the alloy are widely separated in the periodic 
table, the various atomic scattering factors of the individual elements must 
be averaged rather than the atomic number. This program was written for 
copper nickel alloys where the atomic number averaging is a good approxi- 
mation. If the atomic numbers of the constituents of the alloy are not close 
together, the correction to the program is not simple and it may be best to 
skip this correction. The correction is never large for narrow lines and 
usually only makes minor changes in the line shape. The correction can be 
skipped by removing line 40. 
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There are two other line shape corrections that are sotnetiines mentioned in the 
literature. If the path length of the X-rays in the sample depends on theta, differ- 
ential absorption of X-rays takes place as a line is scanned. For the Bragg focusing 
geometry used in a diffractometer, absorption in the specimen attenuates all lines 
uniformly so no absorption correction is needed. A'.x>ther correction is the tempera- 
ture factor which chants line shape because planes with narrower spacing (that is 
greater Bragg angles) are more affected by thermal vibrations than are widely spaced 
planes where the spacing is largie compared to the thermal vibrations. This correction 
is usually the smallest and its calculation is difficult for alloys. It has not yet been 
incorporated into this program. 


137 to These are the atomic scattering factors [9] of pure copper for one-tenth 

143. integer values of sin 9/X. 

145. The atomic scattering factor is calculated for an alloy of copper and nickel 
by averaging the atomic number of the alloy. 

149. The one-tenth integer value of sin 9/ A is calculated for the alloy. This 
value is needed because the table of atomic scattering factors in lines 137 
to 143 gives the factors for discrete values of sin 9/X. If some other 
radiation than copper is used, the factor 1.5417 should be replaced by the 
wavelength of the radiation in angstroms. 

150. This is a linear interpolation between entries in the values in the atomic 

scattering list. is the final correction factor. 

152. If this line is removed, the corrected value of the line will be printed out. 


"pkfd" 

This progfram is located on track 0, file 2 and fits a parabola to three pairs of 
angles and counts that are manually entered. The peak of the parabola is calculated 
and output as the peak value of the line. The three angles must be equally spaced. 
The program is shown in ring 4. 

8 to 13. The I^orentz Polarization Correction is made in exactly the same way that it 
is made in the program "Xfoure '. See the documentation for that program 
for details. 

14 to 35. The values for the atomic scattering factor of both copper and nickel are 
entered for discrete values of sin 0/X. 

36. The atomic scattering factor of the alloy corresponding to the nearest 
discrete value of sin 9/X is calculated. This calculation is slightly more 
accurate than the calculation made in "Xfoure" because the assumption of 
a constant value of the copper to nickel atomic absorption factor ratio is 
not made. 

37, 40. The final values of the etomic scattering factor corrections are calculated 
and 43. using linear interpolation between the discrete values calculated in line 36. 
These factors are for each of the three entered angles. 
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44 to 46. The values of the corrected counts arc printed out so that the operator can 
verify that they are sufficiently close to the peak that the parabola is a 
good approximation. 


”mpkf" 

This program is located on track 0, file 3 and reads pairs of anglM and oounta 
directly from the diffractometer and chooses several sets of three equally spaced 
points to calculate and output the peak from each of these three sets of values. The 
listing is shown in Figure 5. 

1 to 14. The data are read in exactly tire same way that it is read in the program 
"datrd". Angles are stored in the vector A, and counts are stored in the 
vector C. 

15. T is the loop counter which determines the number of sets of points that 
will be used to calculate parabolas. As the program is written, four seta 
will be chosen but any number less than L/2 can be chosen, where L is 
the number of data points entered. 

16. The center point of the parabola is chosen to be which is the largest 
value of the count vector. 

The remainder of the program proceeds exactly the same as *'pkfd”. 


CONCLUSIONS AND RECOMMENDATIONS 


These programs greatly simplify the collection and analysis of data taken from 
the X-ray diffractometer. In addition, a versatile instrument system has been estab- 
lished that an engineer can program to his/her specific needs. 

A number of additional features can be incorporated into the present data 
analysis programs such as: 

1) An automatic decision making program that locates the start and finish of 
the peak. 

2) Calculation of the peak centroid. 

3) A peak location algorithm based on minimization of the asymmetric Fourier 
ooeffleients. 

4) A quick measure of the peak breadth such as calculation of the peak breadth, 
at half height. 
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A 

003 4 '37 

B 

000200 

c 

00*i328 

It 

020000 

A 

00349A 

B 

000200 

f" 

0743:.'? 

li 

020000 

A 

003483 

B 

000200 

c 

004330 

B 

020000 

A 

003508 

B 

000200 

c 

004331 

B 

020000 

A 

003500 

B 

000200 

c 

004332 

B 

020000 

A 

003449 

B 

000200 

c 

004333 

B 

020000 

A 

003399 

B 

000200 

c 

004334 

B 

020000 

A 

003300 

B 

000200 

c 

004335 

B 

020000 

A 

003262 

B 

000200 

c 

004336 

B 

020000 

A 

003141 

B 

000200 

c 

004337 

B 

020000 

A 

003072 

B 

000200 

c 

004338 

B 

020000 


Figure 1. Output from DEC writer. 


fei; 

dip 

''Tills progr am read angles 

1 : 

enp 

"How fiiany data points^", M 

2: 

d 1 m 

ft[M];dirrt CCM] 

3; 

for 

I»1 to t1 

4; 

red 

1 1 , II 

5; 

red 

1 1 , C C I ] 

6 ; 

prt 

” c n ] = " , c c 1 3 

7: 

red 

1 1 , £ 

3: 

red 

11,RE13 

3 ; 

. 0 1 fl [ I ] :: ft C I ] 

10 

: prt 

" ft C I ] = % ft [ 13 

1 1 

; »3rt 

" I = " , I 


12: prt " 

13; next I 

14; files ft t e 5 1 , C t 6 j 1 
15: sprt I, RE*] 

16; spft 
17; end 
-2736 


and counts and stores them on disk. 


Figure 2. Program listing for "dntrd" 
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y: enp "Hom iriiri',' 2 theta 
1 : dim fif C 6 ] } d 1 Hi C #[ & 1 
2i deg 

dim fiCriljdim CCM] 

: #np "Which angle ♦i)e‘’",flf 
trip "Which count file^",Cf 
: Aign 
! *£gri Rf , 1 
8t sread 
9i sread 1,RC»] 

10: f tn t c 3 , b , 5 X , c 7 , / 

11: urt Itf., "2" , 16, "countt" 

12: for 1*1 to M 

13: fmt f 3 . 0, f 6. 2, f 6. 0 

14; I.' r t 1 6 , I , M [ I ] , C t I 3 


15: 

If.: 

17; 

1 I 
1'^; 
20: 

2 1 ; 



34 

35 

.'■t 


39 

40 
4 1 
42 
4 3 

44 

45 
4f 


ne . t I 

d^fj "We ate tioi.i it art wig the background subtract ion 
beep 

enp " Mumbe r i- i.ihere line begins and ends",L,R 
d £ p " P - L fit u s t be € e n " 

C C L 3 J 0 

for I = L t o P 

r [ 1 3-Ci + ( I-L 1 <R-L>*<;C'-C[R3 > >CC 1 3 
If C tl 3 < 0 1 0 > C [ I 3 
f m t f 3 . 0 , f 6 . 0 , f 6 . 2 
gto 27 

i.irt 16, I,CCI3,flC13 
n e X t I 

d £ p " 3 1 ar ♦ i n g L f c o r r e c t i o n " 

dip "Cu radiation and graphite mono, onln." 
prt "LP COPPECTPD PERk IS" 
for I*L to R 

£ 1 n ( H t I 3 2 ' '■ 2c oi R C I 1 2 > ' < 1 ♦ . 75c osfRtI3) '2>>Q 
Q*CL I 3 ICC I 3 

f fii ? f . 0 , f 6 . 2 , , f 8 . 4 

gto 3 3 

ui r t 1 6 , I , H C I 3 , C [ I 3 

prt " 
r I e ;■ t I 
beep 

gsb "RSFCor" 

prt "The integrated intensitM is" 

0 I Q 

f or I = L t o R 
C [ I 3 I.' ; 0 

ne ♦ I 
f d 2 


47 

p'-t 

0 




4 0 

pr t 

" T fi« 1 n» 

e r " a 1 

1 i " , RCR+1 3- 

RtR3 

49 

dip 

" W e ar e 

n o w 1 

ocating the 

peak 

50 

be e p 





51 

e n p 

L h ij O i ► 

t h t e « 

•igual" ipac 

C d 

52 

0 ; T 





5 

pr » 

" C t 0 3 = " . 

C [ 0 3 



54 

prt 

" f C N 3 = " . 

C [H3 



55 

prt 

" i: [ 0 3 * " . 

r C03 




, "degrees two theta" 
«* 


po 1 fit i £ p arit'i 1 ng pe at • 


R-L 


Figure 3. Program listing for "Xfoui’e”. 


e '■ e n " , Q , N , 0 
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56: flCNl-H[Ql>r3 

57: 1 CtO]-l 'CCNl>r2 

58: 1 C[Q]-l'^CCNl>rl 

58 : ft [ Q ) ♦r 3 < < 3*r I +r 2 • ' < 2r 1 +2r 2 > > } P 

66 : f xd 3 

61: prt "The peak is at " , P, "degr ee s two thetsk" 

62: dsp "We are now doing a Rachinger correction" 

63: dsp "this is for a Cu tube only" 

64: 1.54433>r2j 1.54051 in 

65; 3*57. 2958* < <r2-r 1 >t an < ftC L 3 2 > " ( 2r I +r2 > > > ft 

66: prt "Delta theta =",fl 

67: int<2ft '<ft[P+l]-fiCP3)>>D 

68: prt "steps back is",D 

69: dsp "ftssune alpha 1:2 ratio is 2:1" 

70; for I=L+D to R 

71: 57.2958*.00382*t antftC I ] ■ 2 )x l . 54 1 78 ) ft 
72; int (2fl <ftCP + l 3-ft[P] > > >Q 

73; Ct I ]-.5CCI-G!3-.5frc<2ft <ftCP+n-ftEP3>MC[ I-Q-n-CC I-Q]>)Cn ] 

74; if CC I3<0;0>CC 13 

75; next I 

76: for I=R-D to R 

77; 0>C[I3;next I 

78; gto 88 

79: prt "The K alpha 1 line is" 

80; prt "I ftngle Count" 

81: fmt f3.0,f6.2 
82: for I=L to P~D 
83: wrt 16, I ,flC I 3 
84: fxd 5 

85; prt " C C I 3 = " , C C 13 
86: prt " 

87: next I 

83: dsp "How find the K alpha 1 peak" 

89: T>Q 

90; l.'-C[03-l C CN3 >r2 

91: 1. CCQ3-1- CEN3 ; r 1 

92 : ft [ Q 3 +r 3 ( 3r 1 +r 2 > • 2r 1 +2r 2 > i > P 

93; fxd 3 

94: t "The K alpha 1 peak is at ",P, "degrees two theta" 

95: fxd 5 
96; 121H 

97; for 1 = 1 to M- 1 ; fit I 3-P: ftC 1 3 ; next I 

98: for 1 = 1 to P-L+ 1 ; ft C L- 1 + I 3 > ftt 1 3 ; ne t I 

99: ft[13;H 

100: ftC23-ft[ 1 3 >W 

101; P-L+1-D>K 

102: if f re 'K o j=0;l +1 ;k 

103: d 1 111 i' C R ~ L + 1 3 

104; for I=L to R ; C t I 3 ; V 1 1 -L+ 1 3 

105: next I 

106: gto 111 

107: for 1 = 1 to R - L + 1 

108; prt "I=",I 

109: prt "YEI3 = ",Yn3 

110; next I 

111; gsb "FOURE" 


Figure 3. (Continued) 
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112: 

fxd 4 


113; 

prt "Coef ♦' 1 c 1 *r(t s " 


114: 

prt “ft' 0 ■ = " , HI 01 


115: 

for I=l to N 


1 16; 

fmt c 3, f 2. 0, c2, , f 9. 4, , c3, f 2. 0, c2, X, f9 

.4, 

117: 

wrt 16, “fl < “ , I , •■ > = “ ,Ht I 3 Ht0] , "B c “ , I , •* 

>»",B[I3xHC03 

118: 

n«xt I 


119: 

«nd 


128: 

“F0URE”;r4d}crg 2;diM HCOt N] , BCNl | U J; 0> I : <K- 1 >M>D 

121: 

if K 2»int<;K'2>}sfg 2(dsp "ODD NO OF PTS"|r*t 

122: 

If ( I + l > I > >Kjgto -t-e 


123: 

H+ < I - 1 )l4iX j Y[ I 1 > Y; i f I»1 or I«K;l>t1)gt o 

+3 

124: 

If M=4j2>M;gto +2 


125: 

4>M 


126: 

COS' 2;X D)E> >S>F; sin<E> >2>G 


127: 

MYF+HC J3>Ht .T];MVG + BC J] ;B1 Jl 


125:;; 

if ( J + 1 > J ) ; *N; SF -ZG> T } 2F+SG> G; T ) F} gt o - 

1 

129; 

1 > J ; gt o -7 


1 30: 

Yt 1 ]*4Yt2] + Y[|c ]>S;3> I 


1 1 : 

3 4-2Y L I 1 +4Y [ I + 1 1 > 3 ; nnp < 1+-2>I '> = 1 


1 id: 

3W 3D .!H[01; 1> J 


133; 

2 H : n W 3 D .! H [ J 3 ; 2 E [ J 3 H 3 D > B t J 3 ; .1 «i p •: J + 1 > 


134; 

ret- 


1 35 ; 

" flSFCor " ; enp "Nhat atorftic percentage of 

n 1 cl: e 1 is in 

1 36; 

dim F [83 


137; 

29>FC 1 3 


13S; 

25.9>FC23 


139; 

2 1 . 6 ) F [ 3 3 


140; 

1 7 . 9 > F [ 4 3 


141 : 

1 5 . 2 i F [ 5 3 


142; 

1 3 . 3 > F C 6 ] 


143; 

1 1 . 7 : F t 7 3 


144: 

for 1=1 to 7 


145: 

FC I 3*. 01Q*. 0406FC I 3 : Ft I 3 


146; 

next I 


147; 

prt "RTOMIC SCfiTTERING FfiCTOR CORRECTED 

LINE IS" 

143: 

for !=■_ to R 


149; 

1 + 1 n t > 1 0s 1 n H 1 1 3 2 • 1 . 54178 ' > r 4 


150: 

F[r43-<F[r43-F[r4+l 3 >frc< 10sin(,fi[ 1 I- 2). 

1.54178>>H 

151 : 

Ct I 3 H 2>Ct 1 3 


152; 

gto 158 


153; 

fxd 5 


154: 

prt "C [ I 3 = " ,CC I 3 


155: 

prt "H[I3=",ft[I3 


156: 

prt. ■'I = M 


157; 

prt " " 


153: 

ne t I 


159: 

ret 



-» 3€7S 


sampl €?" 


Q 


Fiffure 3. (Concluded) 
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QUigiiiiia. PM& IS 

01 dip "Thii program us«i t h» ♦« point* on * parAbo1« to d*t#r»in# th* p**k“ 

li enp "ang. l",C 

2t tnp "ang. 2",E 

3t »np "ang. 3", I 

4 1 *np "fU)",G 

5i #np "f<2)",H 

6i enp "f <3>" , J 

7: enp "what atomic perctntagt of nickl« i* in th* *anp1*'^",Q 

8 s s 1 n( C t > *» i n( C- 2 > *c os< C ' 2 > #G/ I ♦ . 75c oi<C)*co»<C)>>G 

9s si n(E 2>*»i n<E 2>*co*<>E''2)*H.'< l + . 75 co»(E)*co»<E) ) >H 

101 s i n < I 2 ) » * i n < I 2 > *c o» < I -'2 ) * J ‘ < I ♦ . 75c o* < I ) *c o» < I ) ) > J 

111 since. 2)- 1 .541 78 >M 

12s sinCE 2) 1.54178.>R 

!3s sinC I. 2.>' 1 .54178->S 

14s dim Fl21] 

15s 29>F[1] 

16s 25.9.>FC2] 

17s 21.6>FC3] 

18s 17.9>FC41 
I9s l5.2>Ft5] 

20s 13.3>FC61 
21s 11.7>FC73 
22s 10.2)FC8] 

23s 9.1VFC93 
24i 8.1>F[10] 

25s 23>F[11] 

26$ 25>FC123 
27s 20.75FC13] 

28$ 17.2>FC14] 

29s 14.6fFC15] 

30s 12.7>FC16] 

31s 11.2TFC17] 

32s 9.8>FC183 
33i 8.7>Ftl9] 

34s 7.7>Ft20] 

35s 1 + 1 nt C 1 0*M > 

36s Ftr2]-.01*O*<F[r2]-F[r2*10]>>F[r21 
37s Ftr2 3- <F[r 21-Ftr 2+1 3 M rc C 10M>>K 
38 s 1 + i nt ‘ 1 0F • > r 3 

39s F[r3]-.01Q+' FCrSl -Ftr3+10]> >FCr3] 

40s FCr 33-CFtr 3]-F[r 3+1 3 >f rc C 10R> >L 
41s 1 + 1 nt ( 1 0S • ! t'4 

42: F[r4 3-. 01G)*(FCr4 3-F[f- 4 + 103 > lFtr4 3 

43s F[r4 3- kFtr 43-Ftr 4+1 3 >♦ rc C 10S> >N 

44 s G K 2.>G 

45s H L'2>H 

46s J N 2>J 

47s prt "G*",G 

48s prt 

49s prt 

50s fxd 3 

51s C + C <E-C> -2^ C4H- 3G-J ' C2H-G- J) ) >P 
52s prt "The peak is a» " 

53s prt P 

54s prt "deg. two-theta" 

55s end 
♦ 8346 

Figure 4. Program listing for "pkfd”. 
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M : d 1 hi F I ? 1 ] 

1; «fip "How mano data point 
; «np "Wtiat atoiuic parcent of nickel?”, Q 
duii H[L}idtm C[L] 

4: for 1=1 to L 
red ll,P 
ti red 11,CCI] 
r: prt "C[Il«",Cn3 
St red ll.B 
9: red It, fit I 3 
10: . 0 1 fl C 13 > FI n 3 
It: prt "HCI3»",fltI3 
U: prt "I = ",I 
13: prt " 

14: next 1 

15: for T = 1 to 4 

lb: nia>vCt*3''iH 

IT: for 1 = 1 to L 

IS: if C[I3»H;i3to "LABEL” 

19: ne st I 

J O: "LABEL " : R[ 1 3 >E 

21: ma *C[# 3 >>H 

22: CCI-T31G 

2 5: Ctl+TDJ 

24: flCI-T3)C 

25: ACI+T3>I 

2S: sin>-C 2>*»in',C 2‘*«:os^C 2''*G (l + .r5co*<C>*coi'^C>>)G 

2V; * 1 m E 2 > ♦ ? i n < E 2 ' *c os^ ^ E 2 ' ■>'H < 1 + . 75c os c E > *0 o* >; E > i > H 

2S; s 1 n ' I 2 > 1 n i; I 2 • *c os < I 2 ) ♦ J ' <l + .75cos<I)#cos<l>>>J 

29: sin> C 2>. 1.54178>M 

30: sin<E 2>- 1.54 178 >R 

31: s:nU 2> 1.54 178 >S 

32: 29>Ftl3 

33: 25. 9! Ft 23 

34: 21.6JFC33 

35: 17.9>FC43 

36: 15.2!Ft53 

37: 13.31FC63 

38: 11.7fFC73 

39: 10.2>FC83 

40: 9.1;FC93 

4 1 I 3. 1 >F C 103 

42: 28>Fttl3 

43: 25;FC123 

44: 20.7iFtl33 

45: 17.2>Ftl43 

46; 14.6SFC153 

47: l2.7:Ftl63 

48: 11.2>F[173 

49: 9.S>FC183 

50: 8.7>Ftl93 

51: -.7!FC203 

52: 1 + 1 ft t ■ 1 0 ♦ M > : r 2 

53: F[r2 3-.01bO*<FCr2 3-Ftr2+103 ' >FCr23 
54: F[r23-<F[r2 3-FCr2-M 3 «frc 10M>)K 
55: 1 : nt ^ 1 OR ’ > r 3 

56: Ftr 33-.01O*' Ftr 53 -FCi 5+103i^Ftr33 


Fiffure 5. Program listing for "mpkf" 
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57i Ftr31-<FCr31-FIr3U] >frc(l0l?))L 

58i l+int< ies>>r4 

59i Ftr41-.0lQ*<F[r41-Ftr4*101)>Ftr43 

601 Ftr43-<Ftr4]-Ftr4*U >frc< 10S>)N 

61 1 C'K'£>G 

62i H/L''2>H 

63» J/N'2>J 

64i prt "G»",G 

65i prt 

66i prt 

67i fxd 3 

68i C+<<E-C)/2>< (4H-3G-J) - <2H-G-J)))P 
69 I prt “Th# pe^k Is *t " 

701 prt P 

71i prt "d*g. two-thst*" 

72 1 prt " 

73 I next T 
74: end 
*27620 


Figure 5. (Concluded) 
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